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INTRODUCTION 


Large, flexible orbiting systems have been proposed for possible use in com^ 
munications, electronic orbital-based mail systems, and solar energy collection. The 
size and low weight-to-area ratio of such systems indicate that system flexibility is 
now the main consideration in the dynamics and control problem. For such large, 
flexible systems, both orientation and surface shape control will often be required. 

Figure 1 illustrates a conceptual development plan of a system software capa- 
bility for use in the analysis of the dynamics and control of large space structures 
technology (LSST) systems. This concept can be subdivided into four different 
stages: (1) system dynamics; (2) structural dynamics; (3) application of control 

algorithms; and (4) simulation of environmental disturbances. Modeling the system 
dynamics of such systems in orbit is the most fundamental component. 


SOLAR RADIATION PRESSURE EFFECTS 

The equations for determining the effects of solar radiation pressure on a 
flexible beam are summarized below. 

Forces : 


= -V / 




(t • n) ds (absorbing surface) 


" ^ ^ 2 

Fy = -2h^ / n(x • n) ds (reflecting surface) 


+ e(F^ - F^) 


(surface with reflectivity e) 


Moments : 


/X / 


(absorbing surface) 
s 


My = 2h^ J' n X r(t * n)^ ds (reflecting surface) 


M^ = M^ + e(My - M^) 


(surface with reflectivity e) 
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where 


T =s unit vector in the direction of solar radiation 
/\ 

n = unit vector normal to the surface ds 

h = 4.64 X 10“^ N/m^ 
o 

Results for a flexible beam; 




M„ = 


F,, = 


Y 




-hpa^bQ^Cz^ + Zq)!1 - 2/fljj |a(cos - cos + sin + sin 
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where 

a = sin 6 b^ = cos 0 
o o 

z(x) = flexural deflection 



= nth modal frequency 
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Figure 2 shows the variation of the resultant horizontal and normal force com- 
ponents of a beam with a completely absorbing surface as the solar incidence angle 6 
is varied from 0^ to 90^. Here, 6 represents the angle between the normal to the 
undeflected beam and T. The horizontal and normal force components are measured 
relative to the beam^s undeflected axes. As expected, for small tip deflections of 
the beam, the resultant horizontal absorbing force component becomes zero for inci- 
dence angles of 0^ and 90^, while the normal component has a maximum amplitude at 
zero incidence angle. In figures 2 and 3 the individual effect of each mode is 
illustrated, with the assumed beam tip deflection as indicated in the figures. 

The magnitude of the resultant moments as the solar incidence angle is varied is 
shown in figure 3 for the assumed tip deflection of 0.01£. Large moments can result 
for larger deflections, whereas these moments would be zero for a rigid beam. For 
small pitch angle displacements, the moment due to solar radiation pressure may 
become greater than the moment due to the gravity-gradient forces, as shown in fig- 
ure 4. It is seen that at geosynchronous altitudes, the moment due to solar radia- 
tion may become predominant even for deflections on the order of 0.01£. With the aid 
of moment diagrams such as those in figure 3, it is possible to determine the dis- 
turbance torques due to solar radiation pressure once the number of modes and the 
associated modal deflections are specified for a model. 


MODELING ERRORS - ORBITAL AND GRAVITY- GRADIENT EFFECTS 

One of the principal sources of (disturbance) torques acting on an orbiting 
space structure is the orbital (gyroscopic) and gravity-gradient effects. Such 
effects associated with the orbital (angular) motion do not need to be considered 
when developing a system model for an Earth-based large flexible system. Many 
investigators, however, model the pitch, roll, and yaw modes (rigid body motion) of 
large, flexible orbiting systems as double integrator plants (two poles at the 
origin), and the subsequent control system design is based on these models. It is 
the purpose of this section to evaluate the effects of omitting the orbital and 
gravity-gradient effects when designing shape and orientation control laws for 
flexible systems in orbit. Models of flexible square plates and shallow spherical 
shells in orbit are selected as examples. 

The effects of designing control laws without the orbital and gravity-gradient 
torques included in the system models of square plates and shallow spherical shells 
in low Earth orbit (250 nautical miles) are illustrated in figures 5 to 10. A square 
plate was also considered in a geosynchronous orbit. The analysis was performed by 
first calculating the control law, which was in the form u = Fx, for the case where 
the orbital and gravity-gradient effects are not included in the model. The same 
control law is then applied to the model that includes these effects. 

For figure 4, the control law was selected such that the overall response time 
constant of the system is 2.22 hours (which may be reasonable for a large space 
structure). The shift in the closed loop poles of the plate model due to the pres- 
ence of the orbital and gravity-gradient effects is illustrated in this figure. It 
can be seen that some of the poles move to the right half plane, leading to insta- 
bility and thus emphasizing the importance of including the orbital and gravity- 
gradient effects in the model. The poles due to the rigid body modes are shifted 
considerably, but the flexible modes remain virtually unaffected. This can be 
attributed to the high frequency of the flexible modes. (Note that the orbital and 
gravity-gradient effects are of a relatively low frequency.) This result gives an 
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indication that by designing a more robust (faster response) control system, the 
shift of the actual closed loop poles would be relatively less pronounced. 

This phenomenon can be demonstrated by designing the gain matrix, F, such that 
the desired response time constant is reduced to 460 seconds. The shift in the 
closed loop poles for this case is shown in figure 6. In figure 7 the control forces 
are shown for the second (more robust) control law, where the closed loop response 
with orbital and gravity-gradient effects is degraded but does not become unstable. 

(It should be noted that time has been nondimensionalized with respect to the orbital 
frequency of the 250-nautical-mile low Earth orbit, in order to provide a basis for 
comparison.) The difference in the total control force impulse as applied to the two 
models (a) and (b) is minimal because of the robustness of the controller. The 
slowly varying orbital and gravity-gradient torques have a relatively greater impact 
on the less robust systems, and can even lead to possible instabilities, as was 
illustrated in figure 5. 

As expected, the orbital and gravity-gradient effects are less pronounced in the 
case of a structure in geosynchronous orbit than in the case of a structure in low 
Earth orbit. However, if the control systems are designed with response times com- 
parable to the orbital periods, under the influence of orbital and gravity-gradient 
effects the closed loop systems may become unstable. 

The shift in the closed loop poles of the spherical shell model due to the pres- 
ence of the orbital and gravity-gradient effects is shown in figure 8. One of the 
closed loop poles is moved to the right-hand side of the S-plane, causing instability. 
As compared to the case of the plate, the effect of the orbital and gravity-gradient 
torques on the shell is more pronounced, as the instability due to movement of the 
poles occurs at the relatively fast designed response time constant of 615 seconds 
(compared to 8000 seconds in the case of plate) . When the control is redesigned for 
a response time constant of 400 seconds, the shift in the poles is as shown in fig- 
ure 9. A general shift in the rigid body motion poles is observed, but the system 
remains stable. 

The control forces associated with both models (a) and (b) of figure 9 are com- 
pared in figure 10. A considerable increase in the control effort is observed when 
the model includes the effect of the orbital gyroscope and gravity-gradient torques. 
This may be explained by the fact that the mass distribution of the shell is more 
complex than that of the plate, resulting in relatively greater dynamic coupling when 
the gyroscopic and gravity-gradient effects are included in the shell model. 


THE DEVELOPMENT OF AN ALGORITHM TO EVALUATE COUPLING COEFFICIENTS 

FOR A LARGE FLEXIBLE ANTENNA 

The generic mode equations and the equations of rotational motion of a flexible 
orbiting body contain both coupling terms between the rigid and flexible modes and 
terms due to the coupling within the flexible modes that are assumed to be small and 
thus are usually neglected when a finite element analysis of the d 3 mamics of the 
system is undertaken. In this section a computational algorithm is developed which 
permits the evaluation of the coefficients in these coupling terms in the equations 
of motion as applied to a finite element model of a hoop/column antenna system 
(ref. 1). 
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Using a Newton-Euler approach, one can express the equations of motion of an 
elemental mass of the system, in the frame moving with the body, as: 

*|a^^ + r + 2w + r + a3Xr + 0 )x (Sx r)j^p dv = |f + e + L(q)/p|p dv (1) 

where 

p = mass per unit volume 
e = external forces per unit mass 

q = elastic transverse displacements of the element of volume 
f = force due to the gravity on the unit mass 

L = the linear operator which when applied to q yields the elastic forces 
acting on the element of volume considered 

r = position vector of element dv 

(A) = inertial angular velocity of the body frame 

== acceleration of the center of mass 
cm 


Equations of Rotational Motion 

The equations of rotational motion of the body are obtained by taking the 
moments of all the external, internal, and inertial forces acting on the body; l.e., 
from equation (1) : 


^r X + r + (2o3 x r) + (o) x r) + (o3 x (o) x r)^ p dv = ^ ^ ^ |l(q)/p + f + ^pdv 


( 2 ) 


One can obtain the following form for the equations of rotational motion: 


n=l n=l n=l 


( 3 ) 


where 


* y* jjo X (w X r^) - (r^ • io) (w x p dv 
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OO /• / .. 

X) “ J X q + 2r^ X (u X q) + X (w X q) + q X (u X 9^) 

V 

- (r^ • ui) Qai X q) - (q • (w X rQ>|p dv 


n=l 


oo 

= y* qp dv X (a 


n=l 


oo 

- f ) + 0 )^A A r X dv 

o' Z-rf n n y o 


n=l 


G^ = fr^ X Mr^ p dv 


°° r 

2 “ J (^o ^ ^ p 


n=l 


C = ^ r X e p dv 

V 


r = Tq + q 

M = matrix operator which when applied to r yields gravity-gradient forces 

a = acceleration of the center of mass 
cm 

f^ = force/mass due to gravity at the undeformed center of mass 

= modal shape vector for the nth mode 
60^ = frequency of the nth mode 

= time-dependent modal amplitude function 


Generic Mode Equations 

The generic mode equation is obtained by taking the modal components of all 
internal, external, and inertial forces acting on the body, i.e.. 




• • 


+ 2a)xr + 0)xr + a)x (ci)Xr) pdv 


>] 


= y* • jL(q)/p + f + Jp dv 


( 4 ) 
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The generic mode equation is obtained in the following form: 


\ + “X + + E = (sn + S 8nm + J An 

m=l \ m=l // 


where 


'^n ~ • w X ?o + • CO X (u X p dv 


OO 

= J' • 10 X q + • CO X q + • CO X (oj X q^ p dv 


m=l 


p dv 


( 5 ) 


S Smn = / 

m=l V 



V 


Here is the inertia coupling between rigid body modes and the nth structural 

mode and is the inertia coupling between the mth and nth structural modes. 


Cartesian Components of the Different Coupling Terms 

The expressions for R, G„, \p^, , g , and g in 

^ n mn n mn 

Cartesian components are presented in this section. The following vectors can be 
expressed in their Cartesian component form as: 


= ? i + 5 j 


+ C^k: 


(0 = CO i + 
X 


“y3 + 
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n=l 




V y ^ z 


X 


where i, j, and k are unit vectors along the body principal axes of inertia in 
the undeformed state, and ? , ? , and ? are the coordinates of a point in the 
undeformed state. ^ Y ^ 

With the use of the component forms of the vectors given above, one can expand 
the various vector expressions given in equations (3) and (5) to obtain 


R=[ja) +(J - J)coa)liH-rjd) +(J -1)03 0)^1 

Lxx y yzj Lyy ^ z^zxj-^ 

+ [j 0) 4* (J — J )(j0 o) ”1 k 

L z z y X X yj 


( 6 ) 


Xn) _ 


+ A - (h 

n \ yy zz / x \ : 



h(-)o) 

- H^^> 0 ) 


' zz 

/ X 

yx y 

ZX 



[(n) jj(n)V 

- 

+ 


xy 

yx / y 

\ ZX 

xz 

/ 

z 

1 

+ 0 ) 0 ) ( 

fiCn) ^ 

H(n)' 

y\ xz 

ZX } 

X z\ 

xy 

yx , 


+ - 0 ) 2 ) 

\ z yj\ yz zy j_ 


( 7 ) 




21 


( 8 ) 
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.(n) 


- *n[<»33 - «22>(“yf * “z”') ‘ * “zx ) 


t - i (h'"> + 

x\ yz zy 


+ 2M»_( 

31\ xy yx I 23\ 


+ M., ( + H 


+ - H 


23V“yy 


(n) 

zz 


)] 


) + 0) + (I) 

/ y\ zx xz / z\ xy yx / 


+ w to 
X y 


+ to to + to to (h^“^ + 

\^xy yx / y Z\ yz zy / z. x\ zx xz / 


( 9 ) 


x\ yy zz J y\ zz xx / z\ xx yy / 

tit = zkL + to + (0 

^mn ^%[jx\ yz zy ) y\ zx xz ) z\ xy yx ^ 

+ A fj) + J) + w 

m[”xV yz zy ) y\ zx xz ) z\ xy yx / 

+ to (0 + to to + to to 

+ “x y^xy yx / y z\ yz zy / z xV zx xz J 

- + lM- to^fL^”"^ + - to2(L<”“> + 

x\ yy zz / yV zz xx / z\ xx yy 


5„ - EE m 

a3 


where 


^mn " 




'6 


a3 J ^ot ^3 


a, 3 = X, y, z or 1, 2, 3 
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When a Is X in H 


(n) 

‘oi3 


when a Is 
for 3. 


z In H 


(n) 


a3 


or 

(mn) 

^3 ’ 

the 

a 

is y 

in 

or 

T (®n) 

a3 ’ 

a 


, the corresponding value of a in M „ Is 1. 


,(n) 
*a3 
Ls 3 


or 


(mn) 

^a3 ’ 


In M 


a3* 


a3 

a is 2 in and 

The same reasoning holds 


The expressions for and are obtained by the cyclic permutation 


of X, y, z In the expression for Q 


(n) 


In equation (7) , and the expressions for 


and are obtained by the cyclic permutation of x, y, z in the 

expression for G^ in equation (9) . 

For a discretized model, the expressions for the volume integrals are replaced 
by the following summations: 


“S’ - Z «a>l (fe”*)!”! 


" t Wii (♦S’)* 

i=l 


(a, 3 = x,y,z) 


( 12 ) 


(13) 


where 

k = total number of discrete masses 

i = index identifying a nodal point 

m. = mass concentrated at the ith node 
1 

£ = coordinates of m, in the undeformed state 

i 
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Figure !•- Development of system software for LSST dynamics analysis. 
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(a) Horizontal component . 
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(b) Normal component . 


Figure 2.- Variation of solar force components with incidence angle. 

Totally absorbing surface - free-free beam (length Z = 100 m) . 
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Figure 3 


- Pitch moment due to solar radiation pressure 
surface). Effect of individual modes in the 


(completely absorbing 
system - free-free beam. 



Figure 4.- Moment due to gravity-gradient force as a function of 
pitch angle (100-m rigid beam) . 








(a) Designed closed loop poles. 
Model does not Include orbi- 
tal and gravity-gradient 
effects. 


Closed loop poles resulting 
from the control law of (a) 
when applied to a model which 
Includes orbital and gravity 
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3) 
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3 rigid + 3 flexible modes 

Figure 6.- Shift In closed loop poles due to orbital and gravity-gradient 
torques. Square plate In 250-n.-mi. orbit. Overall designed 
response time constant of the system = 460 sec. 




Forces 


fj, f*3 


Max. Force Astpllcudes 
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f 
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(There is no appreciable difference In forces between 
model with G.G. and orbital effects and model without 
them.) 

ZZ f. (t) - 261.5 N-sec 
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(Mon-Dlmenslonallzed Time) 


Figure 7.- Time history of control forces. Square plate. 




>270 (a) Designed closed loop poles. 

Model doesn't Include orbital 
and gravity-gradient effects. 




(b) Closed loop poles resulting 
from the control law of (a) 
when applied to a model which 
Includes orbital and gravity- 
gradient effects. 
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Figure 8.- Shift in closed loop poles due to orbital and gravity-gradient 
torques. Shallow spherical shell in 250-n.-ml. orbit. Overall 
designed response time constant of the system = 615.0 sec. 
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designed response time constant of the system = 400 sec. 
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PEAK FORCE AMPLITI]DES-C0NTR0L LAW DEVELOPED WITHODT 
ORBITAL AMD GRAVITY' GRADIENT EFFECTS IN THE MODEL AMD 
THEM APPLIED TO A MODEL WITH THEM (Shell in orbit). 
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and G.G. Effects 

With 

6« 6* 

Orbital and 
Effects* 
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Figure 10.- Time history of control forces. Shallow spherical shell. 




